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The running QCD coupling effects for gi singlet are studied. We 
predict that asymptotically gf ~ x~ As , with the singlet intercept As « 
0.85. 

1 Introduction 

The flavour nonsinglet g^ and singlet gf components of the spin structure 
function 171 are presently the objects of intensive theoretical investigations. For 
kinematic region of large x, both g^ s and gf were calculated first in Refs. [1] , [2] . 
Since that, the DGLAP approach has become the standard instrument for the 
theoretical description of 171 . This approach provides good agreement between 
experimental data and theoretical predictions. However, despite this agree- 
ment, it is known that from a theoretical point of view the DGLAP equations, 
being extrapolated into the region of small x, are not supposed to work well be- 
cause they neglect terms ~ (a s ln 2 (l/a;)) n , n > 1 in a systematical way. These 
double-logarithmic (DL) contributions are small when x ~ 1 but become very 
important when x -C 1 and should be accounted for to all orders in a s . Total 
resummation of DL terms was done in Refs. [6] for g^ s and in Ref. [7] for gf. 
These calculations proved that g^ s and gf have the Regge behaviour ~ x~ A 
when x — ► 0. However, the QCD coupling in Refs. [6, 7] is kept fixed. There- 
fore the DL intercepts obtained in [6, 7] are proportional to y/a^ fixed at an 
unknown scale, which makes the results of Refs. [6, 7] to be unclear for practi- 
cal use. Similarly to the DGLAP equations, the DGLAP-like parametrisation 
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a s = a s {Q 2 ) should not be used for gi at small x (see Ref. [8]). Below we 
account for the running a s effects for the description of the small- x behaviour 
of the flavour singlet component of g\. This paper in based on results obtained 
in Ref. [10]. 

2 IREE for the structure function gi 

In order to obtain gi with all DL contributions accounted for at small x, we do 
not use the DGLAP equations. Instead of that, we construct and solve a set 
of infrared evolution equations (IREE), i.e. equations for the evolution with 
respect to the infrared cut-off \i in the transverse momentum space. The cut- 
off /i is introduced as the starting point of the evolution with respect to both 
Q 2 and x. In contrast to the DGLAP equations where only the Q 2 -evolution 
is studied, the IREE are two-dimensional. In order to derive these IREE, 
it is convenient to operate with the spin-dependent invariant amplitude M q 
corresponding to the spin-dependent part of the forward Compton scattering 
amplitude M M! ,(s, Q 2 ) related to g\ as follows: 

gi(s,Q 2 ) = -Z s M q ( Sl Q 2 ). (1) 

7T 

We have used here the standard notations x = Q 2 /2pq and Q 2 = —q 2 > 0, 
with q being the momentum of the off-shell photon and p being the momentum 
of the (nearly) on-shcll quark. The subscript q at M q means that the off- 
shell photon is scattered by the quark. We will need also the spin-dependent 
invariant amplitude, M g (s, Q 2 ) of the other forward Compton scattering where 
the off-shell photon is scattered by a nearly on-shell gluon with momentum p. It 
is convenient to use the Sommerfeld- Watson (SW) transform for the scattering 
amplitudes. However, the SW transform is defined for M q , g : 

— ( S /tl 2 T^)F r (LO,Q 2 ) (2) 

where r = q,g and £(w) = [1 — e~ l7rul ]/2 w ittuj/2 is the negative signature 
factor. 

In these terms, the system of IREE for F r (uj,Q 2 ) can be written down as 
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follows: 



{" + -^) F 9 (^y) = ^2[F 9q {u)F q {u,y) + F gg {uj)F g {w,yj\ 



(3) 



where the anomalous dimensions i 7 ^, with i,k = q,g, correspond to the for- 
ward amplitudes for quark and/or gluon scattering, having used the standard 
DGLAP notations for the subscripts. It is convenient to use normalized ampli- 
tudes Hit : Hh,(uj) — (l/8ir 2 )Fik(oj) , that obey the following equations: 



uH qq — 

WHgg = 



Q<qq ~t~ v qq 



H 



"■99 



qq 

V 99 + H gg 



F^qgHgq, 
- Hg q H q g, 



ujH q 

LUH 



a qg + v qg + H qg (H qq + H gg ), 
a gq + v gq + H gq (H qq + H gg ) . 



(4) 



with 



V ^ = —F){0j),m qq = —, Wlgg 



N 



-2N , m qg = —, Trig q 



NC F . (5) 

Furthermore D(ui) in Eq. (5) accounts for the running QCD effects for Vik, 



1 



1 



(p + rf) 2 + it 2 p + rj 



where r] = \n{p 2 / K 2 QCD ) and b = (33 - 2n/)/127r. 

All Born terms in Eqs. (4) are parametrized as follows: 



where 



with 



"•qq 



A(lo)C f 



2A(lo)N 



"■gg 



7T 



m = - 



T] 2 + TT 2 



dpe 



n f A'(u) 
2ir : 



o {p + r]) 2 +Ti 2 - 
a s A'(u)C F 



i r°° dpe-^p i _ i r 

V~Jo Jp + WibJ Q 



dpe 



P + V 



(6) 



(7) 

(8) 
(9) 

(10) 
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The 7r 2 -terms take place when the arguments of a s are time-like. Finally, 
we arrive at the following expression for the singlet g\: 

/*oo J. . r -, 

i=-{l/xy C + (u,)e n +y + C-(u)e n -y (11) 
-zoo 2m 1 J 



where 



ffgg + #gg ± \ (H qq - H gg ) 2 + AH qg H g 



/2 . (12) 



3 Small- a; asymptotics of g\ 



Eqs. (11,12) give an explicit expression for gf , which accounts for both DL 
contributions and running a s effects. When the limit x — > is considered, one 
can neglect the second term in Eq. (11) and simplify the expression for + . 
Indeed, the behaviour of gf in this limit is driven by the leading singularity 
in Eq. (12). The singularities of f2 + arc related to the branching points of the 
square root. Introducing variables as 

bik = dik + Vik, (13) 

one obtains that the leading singularity of fi + (w) is given by the rightmost 
root, wo, of the equation below: 

^ - A(b qq + b gg )cj 2 + 16(b qq b gg - b qg b gq ) = , (14) 

which predicts that 

gf ~(l/xr(Q 2 /^ 2 r o/2 (15) 

when x — > 0. Eq. (14) can be solved numerically. In our approach, lj depends 
on -q = \n(fi 2 / Aq CD ). The result of the numerical calculation for uio = LOoiv) 
is represented by the curve- 1 in Fig. 1. This curve first grows with rj, achieves 
a maximum where approximately loq — 0.82 and smoothly decreases for large 
rj. In other words we have obtained that the intercept depends strongly on the 
infrared cut-off \x for small values of \i and smoothly thereafter. Quite a similar 
situation was occurring in Rcfs. [9] for the intercept of the non-singlet structure 
function g^ s . We suggest a possible explanation for this effect. The cut-off 
H is defined as the starting point in the description of the perturbative evolu- 
tion. Everything that affects the intercept at scales smaller than /x is attributed 
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Figure 1: Dependence on 77 of the rightmost root of Eq. (14), oj . Curve 2 
corresponds to the case when gluon contributions only are taken into account; 
curve 1 is the result of accounting for both gluon and quark contributions. 



to non-leading effects and/or non-perturbative contributions. Had they been 
accounted for, the intercept would have been ^-independent. Without those 
non-leading/non-perturbative effects taken into account, we then observe an 
important /i-dependence, which becomes weaker for large /i. The maximal per- 
turbative contribution to the intercept As of g\ is obtained from the maximal 
value of uoq. Therefore we estimate the intercept as 

A s = max(co (v)) ~ 0.85 . (16) 

Eq. (16) includes the contributions of both virtual quarks and gluons. These 
contributions have opposite signs and partly cancel each other. It is interesting 
to note that when only virtual gluon contributions are taken into account, this 
purely gluonic intercept, A g is given by the maximum of the curve-2 in Fig. 1, 
which is slightly greater than 1. This value exceeds the unitarity limit, similarly 
to the intercept of the LO BFKL, though in a much softer way. Fortunately, 
by including also the contributions of the virtual quarks the intercept decreases 
down to As of Eq. (16), without violating unitarity. 
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